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Abstract 


Understanding the underlying mechanisms on how solutes modify elastic mod- 
uli of metals is essential in numerous areas spanning solid-state physics to mate- 
rials selection in mechanical design. We perform first-principles calculations to 
study the elastic moduli changes relative to а-Ее for bcc Fe-based solid solutions, 
as solute content increases from 0.4 at.% to 1.85 at.%. Besides the “common ex- 
pectation" that the elastic moduli vary linearly with solute content in these homo- 
geneous solid solutions, nonlinear variations are also observed in Fe-Mn system, 
which is attributed to the change in electronic environment of Fe solvent rather 
than in magnetic moment of solute Mn. At the elastic moduli linear regime, so- 
lutes modify the bulk modulus B and the polycrystalline shear modulus С through 
variations in the valence electron density nws at boundary of Wigner-Seitz cells 
and the bonding strength, respectively. Interestingly, it is found that with increas- 
ing solute content, the change rate of л, Уш (Vm: mole volume) increases expo- 


nentially with the change rate of B, while the change rate of Debye temperature 
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increases linearly with the change rate of С. For 3d transition-metal solutes, the 
change rate of G with solute content increases linearly with electron-to-atom ratio 


е/а from Ti to Mn, but decreases linearly with e/a from Co to Cu. 


Keywords: Fe based solid solution; elastic modulus; solute effect; underlying 


mechanism; density functional theory 


1. Introduction 


In modern industries, body-centered cubic (bcc) Fe based alloys are widely 
used for load bearing, i.e., in strain state, because of low costs and good me- 
chanical properties. For the past several decades, ferritic and martensitic steels 
have been investigated for serving under extreme conditions, such as irradiation, 
corrosion and high temperature. Hereinto, reduced activation ferritic and marten- 
sitic steels have been studied as structural materials in the fusion and the next 
generation fission energy systems [1—5]. Ferritic stainless steels have been con- 
sidered as promising candidates for interconnects in solid oxide fuel cells [6—9]. 
Heat-resistant ferritic and martensitic steels have been investigated as superheater 
materials in fossil fuel power plants [10, 11]. Generally speaking, these variants 
of ferritic and martensitic steels are produced by: i) adding alloying elements for 
solid-solution and the second-phase strengthening, ii) dispersing nanoscale oxide 
particles for improving high-temperature strength and irradiation resistance, and 
iii) surface modifying and coating for mitigating oxidation and corrosion. Efforts 
have been made in the research of strain effects on the microstructural evolu- 
tion, and on the ultimate mechanical properties of ferritic and martensitic steels 
for applications under extreme conditions. It was demonstrated that in metallic 
systems, the stability and mobility of solutes are significantly affected by strain 


[12-15], which might play an important role in radiation-induced segregation. It 
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was also indicated that the applied strain greatly changes the oxidation rate of 
metals [16, 17]. Knowing the molecular mechanisms on how the solutes modify 
the elastic properties of a-Fe can be helpful to predict strain states of ferritic and 
martensitic steels serving as load-bearing components. This can be used to guide 
us for better composition design of this class of steels. 

The effects of solutes on elastic properties in metallic solid solutions have 
been extensively studied [18-20, 23-28]. Miedema et al. suggested an empirical 
relationship between bulk modulus (В) and valence electron density (nws) at the 
boundary of Wigner-Seitz cells for elemental metals, i.e., nws~ VB/Vm, where Vm 
is mole volume [19]. Li and Wu found it applicable to binary intermetallic com- 
pounds and suggested the model: В = пъ У [20]. This relation was verified in 
dilute binary alloys by first-principles calculations [21, 22, 27]. It is reasonable 
that the surface area Aws of Wigner-Seitz cell approximately scales as 2/3 power 
of its volume. From this, it might be deduced an interesting relation: Bown Avie, 
which points out that B is determined by the surface properties of Wigner-Seitz 
cells in metals. Also based on first-principles calculations, Hu er al. reported 
that for W based solid solutions, the shear modulus along [111] direction is pos- 
itively related to the electronic charge density between cations in this direction 
[27]. Many studies revealed the correlations between the change of elastic moduli 
of metals due to alloying and the electronic properties, the volumetric variation, 
the relaxation of bond length, as well as the elastic properties of solutes in pure 
states (see a summary in later Table 5) [21, 22, 27-32]. It might be concluded 
that solutes modify the elastic properties of matrix metals by the combination 
effects of volumetric and chemical bonding changes [28, 31, 32]. Theoretical ev- 
idences indicated that the addition of solutes with higher B in pure states usually 
results in higher B of the metallic solid solutions; but when those solutes with 


lower B in pure states are added, the metallic solid solutions usually have lower B 
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[21, 22, 29]. These studies certainly improve our understanding of solute eftect on 
elastic properties of metallic solid solutions. However, the underlying mechanism 
of how solutes modify elastic moduli of a-Fe is still unclear up to now. 

Early experiments revealed a linear relation of elastic modulus versus solute 
content in bcc Fe based solid solutions [33]. Based on experimental results, some 
studies were done on the correlations between elastic moduli and other properties 
for Fe-base solid solutions [28, 34]. It was found that the rates of change of poly- 
crystalline shear modulus С and Young's modulus E with composition depend 
on the change of lattice parameter with composition and upon the position of the 
solute in the periodic table [34]. Ghosh and Olson decomposed the solute effect 
on G into two contributions: the electronic and the volumetric [28]. They found 
a systematic trend in the electronic contribution as a function of electron-to-atom 
ratio (e/a). Additionally, several theoretical studies on elastic moduli of bcc Fe 
based solid solutions have been done [24, 25, 32, 35-37]. By using all-electron 
exact muffin-tin orbital (EMTO) method in combination with coherent potential 
approximation (CPA), Zhang er al. calculated the single-crystal and polycrys- 
talline elastic moduli of ferromagnetic bcc Fe-X (X = AI, Co, Cr, Mg, Mn, Ni, Rh, 
Si, V) random solid solutions [35-37]. Based on first-principles theory with CPA, 
Khmelevska er al. studied the variations of В with solute content in bec Fe, .Х. 
(X=Si, Ge, Sn; 0 « c « 0.25) disordered solid solutions [24]. Their calcula- 
tions including magnetic properties provide an explanation for the experimental 
non-monotonous variation of B with solute content in Fe-Si system. Fellinger ef 
al. proposed an efficient method for computing solute-induced changes in lattice 
parameters and elastic stiffness coefficients C;; of single crystals using density 
functional theory [32]. However, theoretical studies to elucidate the solute effect 
on elastic moduli of bcc Fe based solid solutions are still needed, to check the 


conventional concepts and to provide new insights. 
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In this paper, first-principles calculations are performed to obtain elastic mod- 
uli for 12 bcc Fe based binary solid solutions with solute content from 0.4 at.% to 
1.85 at.%. To construct these Fe based solid solution systems, we dope а wide va- 
riety of substitutional solutes, including nonmetal Si, simple metal Al, a spectrum 
of 3d transition metals (Ti, V, Cr, Mn, Co, Ni and Cu), and heavy refractory met- 
als (Nb, Mo and W), which have important implications on iron and steel research 
either theoretically or for engineering applications. In this work, the evolution of 
electronic structure and magnetism are carefully examined to seek the underlying 
cause for the variation of elastic moduli versus solute concentration. We investi- 
gate the correlations of the solute effect on elastic moduli, with the solute effects 
on volume, chemical bonding and valence electron density in bcc Fe based solid 


solutions, and with some intrinsic properties of the solutes. 


2. Methodology 


АП first-principles calculations within Density Functional Theory (DFT) are 
performed by using Vienna ab initio simulation package (VASP) [38], and the 
detailed settings can be found in Ref. [39]. Bcc Fe based solid solutions with five 
different solute concentrations are implemented by five periodic supercells that 
each has one solute besides Fe atoms. The dimensions of these supercells and the 
corresponding k-point meshes are listed in Table 1. Evidently, solutes are perfectly 
uniform-distributed in bcc Fe matrix in periodic 3x 3x3, 4x4 x4 апа 5x5x5 
supercells. It should be noted that solutes are not *absolutely" uniform-distributed 
in Ее; X; and Ее Ха systems implemented respectively by periodic 3 x 3 x 4 and 
4 x 4 x 3 supercells, since the number density of solutes in [001] direction is not 
equal to that in [100] and [010] directions for the two systems. For example, in 
the case of Fe;; X, system (3 x 3 x 4 supercell), the number density of solutes 


is 1/4u in [001] direction, while it is 1/3u in [100] and [010] directions. Herein 
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“u” denotes unit cell. As a result, Fe;; X, and Ее Ха systems in equilibrium state 
are actually in body-centered tetragonal (bct) lattice in this work. However, these 
bct lattices are verified to be very close to bcc by the following calculations on 
Ее: Ти and FeosW, systems. After full structure optimizations on the two systems 
by simultaneously relaxing the shape, volume and ions, it is found that the length 
ratio of [001] edge versus [100] (or [010]) edge deviates only slightly from 4/3 
and 3/4, respectively for Fe; Ti; and FeosW, systems. The relative deviations are 
both less than 0.1%. Furthermore, when NS1 (5589) strain is applied on [010] 
and [001] edges of the two systems [39], no matter the [001]/[010] length ratio is 
set to 4/3 (for Еет Ти) and 3/4 (for Еео W1), or equilibrium values, the calculated 
С’ (Сад) shows very small variation, and the relative difference is less than 2.0%. 
These results indicate that the elastic properties of the Fe, X; and Feo5X, systems 
can be approximately regarded in bcc symmetry in this work. Therefore, it is a 
safe approximation that Без, X; and Feog5X, systems are both in bcc lattice. For 
simplicity, the length ratio of [001] edge versus [100] (or [010]) edge is set to be 
4/3 and 3/4 for З x 3 x 4 and 4 x 4 x 3 supercells, respectively. 

In this paper, an arithmetic scheme is employed to extract single-crystal elastic 
moduli of bcc Fe based solid solutions from first-principles calculated stresses 
[39]. Table2 presents some components of strains HS1, HS2, NSI and 5589, 
which are applied on all supercells in our calculations. C' is calculated by the 
equation: С” = (022 — Оп) /12(е22 — €11)]. Other equations for calculating single- 
crystal and polycrystalline elastic moduli can be found in Ref. [39] and references 
therein. From Ref. [40], the elastic Debye temperature O is proportional to a mean 
elastic wave velocity vm: Ө = h/k-(3/(4mQ))'/3-v,,, where h is Planck's constant, k 
is Boltzmann's constant, and О is the atomic volume. In isotropic polycrystals, we 
can consider the relations: pv? = B+4G/3, pv? = С and 3/v;, = 1/v} +2/v}, where 


p is the mass density, and v; and v, are the velocities corresponding to longitudinal 
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and transverse elastic waves, respectively. Since the input values of В and С аге 
at OK, the obtained O is equivalent to that associated with lattice specific heat at 


low temperature limit. 


3. Results 


Experimental measurement of the effect of uniform-distributed solute at dilute 
concentration on elastic moduli of a-Fe might be difficult due to the solute-solute 
and solute-defect clusters [32, 41, 42], and the probable solute nitride as well as 
carbide complexes [43]. However, it is convenient to construct a metallic solid so- 
lution with uniform-distributed solute in molecular modelling [44], so that “pure” 
solute effect on elastic moduli of a-Fe can be theoretically calculated. This makes 
possible to disentangle the contributions on the correlations of the solute effect on 
elastic moduli of a-Fe with the solute effects on other properties of a-Fe, and with 
the intrinsic properties of a solute. In this work, the polycrystalline elastic moduli 
at each solute concentration are subtracted by the counterparts of a-Fe simulated 
with the same group of supercell and k-point mesh, to avoid the systematic error 
due to different k-point mesh. The polycrystalline elastic moduli of a-Fe simu- 
lated with different group of supercells and k-point meshes are listed in Table 3. 
The calculated differences, AB, AE and AG of Fe-X (X- 12 solutes) systems, rep- 
resenting the changes of elastic properties owing to solutes, are presented on the 
left panel in Figs. 1 and 2. As it can be seen, for 11 Fe-X (X-Ti, V, Nb, Cr, Mo, 
W, Co, Ni, Cu, Al and Si) systems, these elastic moduli decrease linearly in dif- 
ferent slopes, when solute content increases from 0.4 at.% to 1.85 at.%. As shown 
in Fig. 1(g), for Fe-Mn system, AB increases abruptly and drastically at 1.39 at. 96 
Mn, while AG and AE increase linearly with Mn content. 

The calculated lattice constant changes Aa (Aa = a — ao) owing to solutes, 


are presented on the right panel in Figs. 1 and 2 for all systems. Aa values are 
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calculated by a of Fe based solid solutions subtracting ay of a-Fe simulated with 
the same group of supercell and k-point mesh (see Table 3), to avoid the influence 
of different k-point mesh. As we can see, for 10 Fe-X (X=Ti, V, Nb, Cr, Мо, W, 
Co, Ni, Cu and Al) systems, Aa increases linearly in different slopes when solute 
content c increases from 0.4 at.% to 1.85 at.%. However, ^a of Fe-Mn and Fe-Si 
systems show very small change as c increases. It should be noticed that Aa of 
Fe-Mn system shows an appreciable deviation from the fitting line, at 1.39 at. 96 
Mn. Additionally, the change of Debye temperature, AO, of these Fe-X (X- 12 
solutes) systems, are calculated in a similar way. As shown on the right panel in 
Figs. 1 and 2, AO of Fe-Mn system increases linearly with Mn content, while AO 


of all other systems decreases linearly with solute content. 


4. Discussion 


4.1. Electronic origin of the nonlinear variation of B with Mn content 


It is long established that elastic moduli of bcc Fe based solid solutions with 
uniformly or randomly distributed solute are linear functions of solute content 
within solubility limit [28, 32-34]. This model is well confirmed by our calculated 
AG, AE and AB of 11 dilute bcc Fe-X (X= Ti, V, Nb, Cr, Mo, У, Co, Ni, Cu, AI 
and Si) systems. However, it is not followed by the AB of Fe-Mn system. This 
linear behavior can be understood for atomic configurations with nonexistence of 
solute-clusters, impurities and defects. However, the condition from the electronic 
structure perspective has not been discussed so far. Noticeably, the nonlinear 
jumps of AB and Aa in Fe-Mn system are synchronous, as shown in Figs. 1 (g and 
g’). In the following part of this subsection, we deduce the necessary condition 
in electronic structure for the linear dependence of lattice constant a on solute 
content c in metallic solid solutions, and demonstrate that it is also correct for the 


linear dependence of bulk modulus B on c. 
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In a metallic solid solution M;_.S., the atomic volume Q, defined as the mean 


volume per atom, usually varies in a linear fashion with c: 
О = (1 —ООм + сос. (1) 


Herein, Ом is the atomic volume of matrix metal M, and O% is the effective 
atomic volume of solute S. For some metallic solid solutions, this linear relation 
is only valid below a maximum solute content cmax, since the constant parameter 
O% changes when c > cmax. According to King’s explanation, as solute content 
c increases beyond cmax, sufficient amount of solutes alter the electronic environ- 
ment of the solvent, leading to the change of the effective solute volume Qe [50]. 
This implies that the invariance of the electronic environment of the solvent is a 
necessary condition for the linear dependence of © on c. In а dilute bcc Fe based 


solid solution, the lattice constant a varies linearly with solute content c on the 


condition that c is below a certain limit стах: 
а = ад + сл  CE€Cmax. (2) 


where ад is the lattice constant of a-Fe апа r, is the change rate. Obviously, for 
a dilute metallic solid solution, cr, is very small compared with ag. If one takes 
the third power on both sides of Eq. (2), and omits those terms of higher than one 


power of cr, on the right side, there is 
а -(1— с)аё + c(ag + Зт.)ад. (3) 


After both sides of Ед. (3) are divided by two, it is identical with Ед. (1). We 
obtain: О = а?/2, Ом = аў/2, and Qc = (ao + 3ra)as/2. When c increases 
beyond стах, Г. changes so that a deviates from the linear trend and the linear 
relation of €2 versus c also breaks down according to Eqs. (2) and (3). Therefore, 


the nonlinear variation of a with c has the same electronic origin as that of 0 
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in dilute Fe based solid solutions, and the linear dependence of a on c implies а 
necessary condition in electronic structure as that of О. 

Towards understanding the electronic origin of nonlinear variation of bulk 
modulus В with c in dilute bcc Fe;_.-Mn, system, we examine the evolution of 
differential charge density рат in (110) plane with c in 8 Fe-X (X=Ti, Cr, М, Mn, 
diff ; 


Co, Ni, Cu and S1) systems. Herein, o 15 defined as the difference between the 


total charge density of a system and the superposition of the isolated atomic charge 
density placed at atomic sites [45]. One can observe that the ptf around solute 
Mn changes when its content increases from 1.04 at. % to 1.39 аі. %, as shown 
in Fig. 3. For example, the isodensity olive ring transforms from cross-like shape 
into rhombus-like shape. This clearly shows the electronic environment change 
around Mn, which leads to the nonlinear variations of AB and of Aa with Mn con- 
tent. On the other hand, we could not observe any variation in nmi pictures of 
(110) plane in other systems as solute content c increases, whose AB and Aa are 
both in linear dependence on c. These results confirm the above deduction that in 
homogeneous metallic solid solutions, the invariance of electronic environment of 
solvent is a necessary condition for the linear dependences of B and of a on solute 
content c. Specifically, once sufficient solutes alter the electronic environment 
of bcc Fe solvent, the “sensitive” bulk modulus varies nonlinearly at this solute 
content. 

Now, let us turn our discussion on the perspective of the evolution of mag- 
netism in these bcc Fe based solid solutions as solute content c increases from 
0.4at.% to 1.85at.%. Herein, the magnetic moments of ions are calculated in 
equilibrium volume, by using the standard Wigner-Seitz radii in the VASP pseu- 
dopotential database. The magnetic moments of solutes and Fe ions do not change 


significantly with c in 7 Fe-X (X-Ti, Cr, W, Co, Ni, Cu and S1) systems, except 


in Fe-Mn system. Figure 4(a) displays that Fe ions in the nearest two and farther 
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neighbor shells (Inn, 2nn and farther) around Mn exhibit very small variation of 
magnetic moment и, within this range of Mn content; furthermore, the и values 
of the Inn and 2nn Fe ions are slightly smaller than that of farther Fe ions. These 
results are in good consistence with the experimental conclusion that the magnetic 
moments of neighboring Fe ions are almost not affected by solute Mn [46]. As 
shown in Fig. 4(b), the д of Mn increases monotonously from —1.84 ив to 0.36 
ив, as Mn content increases from 0.4 at.% to 1.85 at.%. The magnetized localized 
state of Mn in ferromagnetic iron was firstly deduced by Jaccarino et al. from 
experiments [47, 48]. The calculated и of Mn is 0.11 ив at 1.39at.% Mn, which 
agrees well with the experimental result, nearly zero, in a bcc Fe-Mn solid solu- 
tion with 1.5 at.% Mn [46-49]. Our calculated д of Mn is —1.56 ив in Ее Ми 
system (0.78 at.% Mn) with lattice constant of 2.8323 А. This result is very close 
to the EMTO result, —1.49 ив, but in a large difference with the PAW-PW91 one, 
0.9 up, in Ее! Мп) system reported by Olsson et al. [41]. We specially set the 
lattice constant to 2.83 À for Ее Ма; system, and the calculated u of Mn, 0.56 
Ив, 15 roughly consistent with the result of Olsson et al. Nonetheless, this volume 
effect on magnetism of Fe-Mn system might need further research. The Aa vari- 
ation with Mn content c іп Fe-Mn system might be attributed to the и change of 
Mn. For example, when c changes from 0.4 at.% to 0.78 at.%, ^a shows a "large" 
increase, by 0.0011 А (see Fig. 4(c)), which might be owing to stronger magnetic 
repulsion among Mn ions due to denser population. Furthermore, it is owing to 
the smallest magnetic repulsion between Mn and other cations, that Fe-Mn system 
has the minimum volume at 1.39 at.% Mn. However, this magneto-volume effect 
may not fully take in account the observed elastic moduli changes in this dilute 
Fe-Mn system, in view of the predicted hydrostatic-strain effect on elastic moduli 
of a-Fe in Ref. [31]. For example, from 1.04 at.76 Mn to 1.39at.% Mn, 1) the 


elastic moduli do not increase uniformly, e.g., 2.20 GPa of AC44 versus 17.18 GPa 
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of АСи, as shown in Fig.4(d), which do not indicate a hydrostatic-strain origin; 
2) the B increment of 17.67 GPa can not be induced by the very small change of 
lattice constant, —0.00185 А, equivalent to hydrostatic-strain decrease of -0.002. 
Therefore, the nonlinear variation of B, Си and Со can only be attributed to the 


change of electronic structure in the dilute bcc Fe-Mn system. 


4.2. Exponential relation between д(п Vj.) /дс and B, 


King defined an atomic volume size factor for metallic solid solutions [50], 


Qy =. (4) 


Herein, с denotes solute content; Q and OM are atomic volumes of metallic solid 
solution and of solvent matrix, respectively. О, expresses the relative change 
rate of atomic volume with solute content c, and characterizes the solute effect on 
atomic volume of the solvent [41, 50]. Similarly to О, у, an elastic modulus factor 


can be defined for bcc Fe based solid solutions: 


1 дк 
Кре дс i 


(5) 


Кор = 


Herein, Kfe and к represent elastic moduli of a-Fe and of bcc Fe based solid 
solution, respectively. Ку expresses the relative change rate of elastic modulus 
with solute content c, and characterizes solute effect on an elastic modulus of a- 
Fe. We will discuss correlations between the solute effect on an elastic modulus of 
a-Fe and other quantities via Kmr in the following part. Since Debye temperature 
© reflects bonding strength of a crystal [51, 52], its “modulus” factor, бъ, is 
proposed to characterize the solute effect on bonding strength of a-Fe. 
Interestingly, let us consider the correlations between the effect of solute on B 
and some other quantities. Figure 5(a) indicates no correlation between В,„ and 
Өм in these dilute bcc Fe based solid solutions. As shown in Fig. 5(b), there is 


a very weak positive correlation between Ви; and Ог, which is inconsistent with 
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the concept that volume expanding leads to decreasing elastic modulus [31, 32, 
39]. However, it should be noted that there are good negative correlation between 
B and equilibrium atomic volume in face-centered cubic (fcc) Al based, bcc W 
based and hexagonal close-packed (hcp) Mg based solid solutions [22, 27, 29]. 
Our results show that solute effects on bonding strength and on volume do not 
contribute much in modifying B of a-Fe. Figure 5(c) shows that the calculated 
Bmp are in a rather good positive correlation with experimental В of solutes in 
pure states [53], which agrees with the reported positive correlation between the 
calculated B of fcc Ni based, fcc Al based, and hcp Mg based solid solutions, and 
experimental B values of pure solutes [21, 22, 29]. 

In order to explain what is behind the positive correlation between B of pure 
solutes and B,,;, we calculate the valence electron density nws in bcc Fe based 
solid solutions. We firstly calculate B with PAW-PW91 functional for all elemen- 
tal solids except for Mn [54], since its ground state, the noncollinear antiferromag- 
netic a-Mn, needs complicated initial settings and time-consuming computations 
to converge into the correct magnetic order [56]. Table4 lists Ше PAW-PW91 
B values of pure solutes where that of Mn is from Ref. [56] also calculated with 
PAW-PW91 functional [41, 58]. Figure 6(a) displays a good positive correlation 
between these theoretical B values of pure solutes and Ви; of Fe based solid solu- 
tions. From previous studies in Refs. [19—22, 27], it is reasonable that В~п Vin 


is applicable to bcc Fe;_.X, solid solutions. Similarly to Ref. [27], it is defined 
nws = ((1 — с)пв Уве + cnx Vx)/ Уп. (6) 


Herein, nx = VBx/Vx with Bx, Vx and nx in units of GPa, 10-5m?/mole and 
e/au?, respectively [20]. Уве, Vx, пв. and nx can be calculated from theoretical 
results of lattice constants and B listed in Tables3 and 4. If one takes partial 


derivative on both sides of the expression В-па Vm with respect to solute content 
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с, there is 


Въ, Вве-д( пък У„)/дс. (7) 


As shown in Fig. 6(b), the theoretical data of д(пух Vin)/Oc (divided by Br.) except 


that of Fe-Nb system can be perfectly fitted to the equation: 
долу Уп)/дс = руехр(рэВ» ў) + ps. (8) 


The adjusted В? is 0.887, and the fitted values of ру, ро and p; are 2.1, 0.268 and 
— 1.241, respectively. Herein, for better display the exponential trend, the data of 
Fe-X (X-Re, Os and Ir) systems are included in Fig. 6(b). As shown in Fig. 7(a), 
the first-principles AB of Fe-X (X-Re, Os and Ir) systems all have gentler slope 
than that of Fe-W system, indicating larger B, than that of Fe-W system. Our 
results reveal an exponential, rather than linear, correlation between alnis Vin)/Oc 
and Ви; in bcc Fe based solid solutions. Noticeably, nws is quite different from 
the concept of free electron density “free, and the evidence is as following. Our 
calculated nws decreases with solute content с in all Fe-X systems. However, the 
first-principles Fermi energy er increases with с in 10 Fe-X (Х=АІ, Cr, Ir, Mn, 
Mo, Os, Re, Si, У and W) systems. This means that пес. also increases with с, 
since £p scales as 2/3 power of пас, according to free electron gas model. Hence, 
nws and ne, Show opposite trends with solute conte с in these systems. 

Let us then, discuss the mechanism of how solutes modify B of a-Fe. As 
common knowledge, B characterizes the compressibility of solids. For metal- 
lic systems, B is mainly determined by the kinetic energy of free electrons and 
the overlapping between bound-electron clouds around neighboring cations which 
contribute repulsion forces [44, 59]. Actually, the repulsion between overlapped 
bound-electron clouds can be attributed to increased kinetic energy of electrons in 
the overlapping area due to increased density. According to Thomas-Fermi statis- 


tics, the kinetic energy of electron cloud scales as 2/3 power of its density [59]. 
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Hence, there is an inverse correlation between В and volume of metals. The gen- 
eral applicability of the empirical relation nws- VB/V, in elemental metals and in 
binary alloys [19-22, 27], illustrates that the so-called valence electron density at 
boundary of Wigner-Seitz cells, nws, perfectly represents the combination effects 
of free electron density and the overlapping between bound-electron clouds in a 
phenomenological way. From this point, it is convenient to discuss the electronic 
origin of solute effect on B of a-Fe Бу nws. 

Figure 8(a) presents the nws on a Wigner-Seitz cell in a-Fe. Solute contributes 
free and bound electrons into the metallic solid solution, so that it directly changes 
nws; furthermore, solute can modulate nws by changing the volume. Figure 7(b) 
presents the first-principles Aa of Fe-X (X=Re, Os, Ir and W) systems versus so- 
lute content, indicating that their volumes are much expanded by adding solutes. 
Our calculated О, ; of solutes can be sorted into two classes. The first class in- 
cludes solutes Ir, Mo, Nb, Os, Re and МУ, whose ©, аге fairly large, e.g., 0.628, 
0.545 and 0.534 respectively for Ir, Os and Re. The second class contains all the 
other solutes, whose О, у is much smaller than those of the first class, as shown 
in Fig. 5(b); for example, the greatest two 2; values among the second class, аге 
0.37 and 0.296, respectively belonging to Ti and Ni. Solutes in the first class have 
much bigger volume-expanding effect than those in the second clase, leading to 
greater decrease of B. This explains the left shift of the half-filled symbols (rep- 
resenting solutes in the first class) from the short-dash line fitting to those filled 
symbols (representing solutes in the second class), and the exponential relation 
between ду Vm)/Oc and Ви; in Fig. 6(b). Hence, solutes modify B of a-Fe pri- 
marily via the valence electron density, nws, at boundary of Wigner-Seitz cells, 
including by modulating nws via volume. Table 5 lists the verified correlations 
between B of metallic solid solutions and other properties in previous studies and 


this work. Obviously, there exists similar correlation between nws and В in hcp 
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Mg-X (X=solutes) solid solution as those in other systems owing to the correla- 
tion between В and Вх. It is reasonable that good linear relation nws~ VB/Vin 
is unavailable in these metallic solid solutions because of the volume effect of 


solutes, as in our case of bcc Fe-X system. 


4.3. Linear relation between Ф. and Ging 


Here we turn to discuss the correlations between solute effect on polycrys- 
talline shear modulus С of a-Fe and other quantities. As shown in Fig. 9(a), On, 
increases in a perfect linear fashion with Су except for Fe-W system. The calcu- 
lated Ө, of Fe-W system is much smaller for matching Су, which is attributed 
to the much heavier mass of W ion. However, if the mass of W was replaced by 
that of Fe in calculating © of Fe-W system, the resulted ОФ; is right located on 
the fitting line, as shown by the open diamond symbol in Fig. 9(a). О, exhibits 
almost no correlation with Соу, as shown in Fig. 9(b). These results indicate that 
the solute effect on bonding strength, rather than on volume, plays a dominant 
role in modifying G of a-Fe. This can be well understood from the "spring" rep- 
resentation of metallic bonds in a-Fe as shown in Fig. 8(b): all these bonds resist 
shear strain in a consistent manner. We also observe similar linear correlation be- 
tween Ө, and E, as that between ©,,¢ and С, ў, and no correlation between €; 
and Ему. There is a weak positive correlation between G, and experimental С of 
pure solutes [53], as shown in Fig. 9(с). Ems also show a weak positive correlation 
with experimental Е of pure solutes. 

Hu et al. found that the shear modulus in [111] direction of bec W based 
solid solutions, Су, is in positive correlation with Ше W-X (X=solutes) bonding 
strength in [111] direction [27]. They reported that С of W-Mo solid solution is 
greater than that of other W-X (X-Ti, Pd) systems, which is corresponding to the 
higher charge density between adjacent W and Mo atoms in [111] direction than 


that between W and Ti (or Pb) atoms. They also reported that the theoretical G 
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values of bcc W based solid solutions, exhibit very weak correlation to equilibrium 
volume [27]. Their results are consistent with our calculations of Fe based solid 
solutions. As far as our reach, there has been no more report on the correlation 
between G and bonding strength in metallic solid solutions, as shown in Table 5. 
Nonetheless, the perfect linear relationship between G,,; and Ө, in bcc Fe-X 
system and the results in bcc W-X system [27] suggest that solutes modify G 
of bcc metals mainly via variation of bonding strength. Additionally, previous 
first-principles calculations indicate G of fcc Ni based solid solutions are in good 


negative correlation to equilibrium volume, which may need further study [21]. 


4.4. Correlation between Ст; and e/a 


From above subsection, it is concluded that solutes modify С of а-Ее mainly 
by modulating electronic interactions between neighboring cations. Herein, we 
further discuss the case of 3d transition-metal solutes. Figure 10 plots Си; versus 
electron-to-atom ratio e/a of 3d transition metals, which provides instructive in- 
formation on how the electronic interactions between Fe and 3d solutes influence 
G of a-Fe. As it can be seen, Су approximately exhibits a linear increasing trend 
with e/a for early 3d metals up to Mn, but a linear decreasing trend for later ones 
from Co, so that the fitting lines present in a volcano's shape as a whole. This indi- 
cates that transition-metal solutes with similar number of valence electrons to iron 
form stronger bonding to bcc Fe solvent, leading to greater G. Hu et al. reported 
similar results for bcc W based solid solutions: an approximately linear increas- 
ing trend of G with the increasing number of valence electrons of early transition 
metals lying to the left of W in periodic table, but an approximately linear de- 
creasing trend for later ones lying to the right of W [27]. Nonetheless, it should be 
noted that the increasing trend of Су versus e/a from Ti to Mn is opposite to the 
linear decreasing trend of дС /дс versus e/a for Fe-X (X=3d transition-metal so- 


lutes) systems in Ref. [28]. This discrepancy might be due to structural difference 
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between the bcc Fe based solid solutions with uniform-distributed solutes in our 
modelling, and real materials, in which solutes almost always exist in short-range 


order [42]. 


5. Summary 


Based on first-principles calculations within the density functional theory, the 
changes of elastic moduli (AB, AG and ЛЕ) owing to solutes, are calculated for 12 
bcc Fe based homogeneous solid solutions within solute content c from 0.4 at.% to 
1.85 at.%. Results show that these elastic moduli decrease linearly with c in 11 Fe- 
X (X=Ti, V, Nb, Cr, Mo, W, Co, Ni, Cu, Al and Si) systems. For Fe-Mn system, 
AB show a nonlinear increase at 1.39 at. % Mn, while AG and AE increase linearly 
with Mn content. It is found that the nonlinear variation of AB in Fe-Mn system 
originates from the change in electronic environment around solute Mn rather than 
in magnetic moment of solute Mn. It is concluded that the elastic moduli of bcc 
Fe based solid solutions are in linear solute-content dependence on the condition 
of unchanged electronic environment of bcc Fe solvent, and vary nonlinearly once 
sufficient solutes change the electronic environment of the solvent. 

In this paper, an elastic modulus factor ки; is defined to characterize solute 
effect on the elastic modulus к of a-Fe. The “modulus” factor of Debye tempera- 
ture, Ө, is proposed to characterize solute effect on bonding strength of a-Fe. It 
is concluded that solutes modify B of a-Fe primarily via the valence electron den- 
sity, nws, at boundary of Wigner-Seitz cells, and modify С via bonding strength. 
It is found that the theoretical Os Vin)/Oc increases exponentially with the bulk 
modulus factor Ву, and ду increases linearly with Су. For 3d transition-metal 
solutes, Gmr increases linearly with electron-to-atom ratio e/a of solutes from Ti 
to Mn, but decreases linearly with e/a of solutes from Co to Cu, indicating that by 


adding solutes having similar valence electrons to iron leads to greater G of bcc 
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Fe based solid solutions. These results contribute new insights on the underlying 


mechanisms of how solutes modify elastic moduli of a-Fe. 
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Figure Captions: 

Fig. 1. (Color online) Changes of bulk modulus, polycrystalline shear modu- 
Ius and Young's modulus, AB, AG and AE (left panel), and lattice constant change 
Ла, Debye temperature change ЛО (right panel) of bcc Fe-X (X= early transition 


metals) sold solutions; the abscissa 1s solute content. Lines are linear fits. 


Fig. 2. (Color online) Changes of bulk modulus, polycrystalline shear modulus 
and Young's modulus, AB, AG and AE (left panel), and lattice constant change 
Ла, Debye temperature change AO (right panel) of bcc Fe-X (X=later transition 
metals, Al, and Si) systems; the abscissa is solute content; all the scales are the 


same as those in Fig. 1. Lines are linear fits. 


Fig. 3. (Color online) Differential charge density on the (110) plane in bcc Fe-Mn 
solid solution with Mn content 1.04 at. % (left panel) and 1.39 at. 90 (right panel). 
In both pictures Mn atoms are placed at the center. The unit of charge density is 


е/АЗ. 


Fig. 4. (Color online) (а) Averaged magnetic moments of the first and second 
nearest neighboring (Inn and 2nn) shells of Fe atoms around Mn, and farther Fe 
atoms, and (b) magnetic moment of solute Mn in bcc Fe-Mn solid solution at 
different Mn contents; (c) lattice constant change Aa, and (d) changes of single- 
crystal elastic moduli АСи, АС and AC and bulk modulus change AB of bcc 


Fe-Mn solid solution at different Mn contents. 
Fig. 5. (Color online) (a) “Modulus” factor of Debye temperature ду, (b) volume 
size factor О, у, and (c) experimental bulk modulus В of pure solutes versus bulk 


modulus factor Ви; of bcc Fe-X (X= 12 solutes) systems. Ви; of Fe-Mn system 
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is calculated in the range from 0.4 to 1.04 at.% Mn. 


Fig. 6. (Color online) First-principles results of (a) bulk modulus B of pure solutes 
and (b) пус Vin)/Oc (see the text for a description) versus bulk modulus factor 
Вт; of bcc Fe-X (X= 15 solutes) systems. Ви; of Fe-Mn system is calculated in 
the range from 0.4 to 1.04 at.76 Mn. Solid curve is exponential fit to filled and 
half-filled symbols; short-dash line is linear fit to filled symbols. 


Fig. 7. (Color online) (a) lattice constant change Aa and (b) bulk modulus change 


AB versus solute content for bcc Fe-X(X= Ir, Os, Re and W) solid solutions. 


Fig. 8. (Color online) (a) Charge density on a Wigner-Seitz cell in a bcc cell 
of a-Fe. The lower two hexagons are darkened to improve stereoscopic impres- 
sion. The dash lines are body diagonals. The unit of charge density 15 е/АЗ. (b) 


"Spring" representation of metallic bonding in a bcc cell of a-Fe. 


Fig. 9. (Color online) (a) “Modulus” factor of Debye temperature Ө, (b) vol- 
ume size factor О, ў, and (c) experimental polycrystalline shear modulus С of pure 
solutes versus the polycrystalline shear modulus factor Су of bcc Fe-X (X= 12 
solutes) systems. In panel (a), the vertical coordinate of the open diamond symbol 
marked Бу W* is the ду value of Fe-W system calculated by replacing the mass 


of W with that of Fe; the line is a linear fit to the data except those of Fe-W system. 
Fig. 10. Polycrystalline shear modulus factor Си; versus electron-to-atom ratio 


е/а of За transition-metal solutes. Lines are linear fits. The electronic configura- 


tions of 3d transition metals in ground state are shown. 
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Table 1: Atomic percent concentrations of solutes (at.%), supercells (in multiples 1, m and n of a 
two-atom bcc Fe cell respectively along [100], [010] and [001] directions) and k-point meshes of 


five model systems. Herein X denotes solute; /, т and n are integers. 


system at.06 Х supercell k-point mesh 


ЕелХ, 0.40 5x5x5 4x44 
Fe;i5;X, 0.78 4х4х4 54555 
FeysX, 1.04 4х4х3 557 
Ее Х, 1.39 3x3x4 7*7*5 
Fes3X1 1.85 3x3x3 7*7*7 


Table 2: Components (%) of HS1, HS2, NS1 and SS89 strains. 


HS1 (ен!) HS2 (ен) №1 (ен, е», ез) SS89 (е, €33, €12) 
0.15 —0.15 —0.99, 1.0, 0 0.004, 0, 0.873 
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Table 3: Lattice constant ag (Unit: А), polycrystalline elastic moduli B, E and С (Unit: GPa) and 


Debye temperature © (Unit: К) of a-Fe calculated by different group of supercells and k-point 


meshes. 
system do B E G о supercell k-point mesh 
Feo; 2.8303 191.50 214.54 81.68 468.57 5х5х5 4»4»4 
Fes 2.8305 192.00 214.61 81.68 468.61  4x4x4 5%5*5 
Ее» 2.8300 193.50 215.76 82.09 469.78 4х4х3 5s5x7 
Fen 2.8297 195.00 216.15 8271 471.53  3x3x4 TXT*5 
Fes, 2.8292 195.17 218.68 83.26 472.99  3x3x3 771 
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Table 4: Lattice constant a (and c for tet and hcp structures; Unit: À) and bulk modulus B (Unit: 
GPa) of elemental crystals calculated in this paper [54]. Herein tet and dia are abbreviated for 
tetragonal and diamond structures, respectively. AFM and FM are abbreviated for antiferromag- 
netic and ferromagnetic phases, respectively. Experimental (exp.) values are also listed for com- 


parison. 


elemental crystal а, c (PW91) a,c(exp)?* B(PW91) B (exp.)> 


V bcc 2.9766 3.02 188.33 158 
Nb bcc 3.3202 3,30 182.56 170 
Cr bcc AFM 2.8482 2.88 196.67 160 
Mo bcc 3.1703 3.15 249.67 261 
W bcc 3.1740 3.16 312.0 311 
Ir fec 3.8789 3.84 345.67 371 
Ni fcc FM 3.5229 3.52 192.56 177 
Cu fcc 3.6285 3.61 133.11 138 
A] fcc 4.0465 4.05 74.78 75,2 
Si dia 5.4665 5.43 89.11 98 

Mn tet AFM 8.669, 8.668° 8.877, 8.8739 1889 92.6 
Ті hep 2.9195, 4.6195 2.95, 4.69 122.56 108 
Re hcp 2.7774, 4.4832 2.76, 4.46 371.17 334 
Os hcp 2.7588, 4.3512 274,433 401.83 373 
Со hcp FM 2.4862, 4.0169 2.51, 4.07 207.67 182 


аве [55]  PRef.[20]  CRef.[56] “Ref. [57] 
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Table 5: Correlations of В and С of metallic solid solutions with other properties. nws is valence 
electron density at boundary of Wigner-Seitz cells and V is volume of metallic solid solutions. BS 
is the acronym for Bonding Strength. Bx and Gx are bulk and polycrystalline shear moduli of 
solute X, respectively. y and x denote that correlation exists and do not exist, respectively; v v 


denotes quantitative relationship. 


nws V Ву BS Gx 
- ¥ | G(fee Ni-X) [21] 
— | G(bcc W-X) [27] N 
G(bcc Fe-X) [28] 

G(hcp Mg-X) [29] 


G(bcc Fe-X) [ours] vv x x 


B(fcc Ni-X) [21] y 
B(fcc Al-X) [22] y 
B(bcc W-X) [27] y 
B(hcp Mg-X) [29] - 
B(bcc Fe-X)[ours] УМ 
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